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SYMBOLS 


C, 

viscous  damping  constant  between  m, 

V 

scaled  velocity  y/c.  for  one-degree-of- 

and  m.  i 

f  reedom  system 

F, 

forc  ing  function  on  m, 

V' 

scaled  velocity  at  t  =  f.  for  one-degree- 
of -freedom-system 

fto 

nonlinear  spring  force  between  mass  and 
foundation;  kx  +  fix 1 

V, 

scaled  velocity  y./w. 

A. 

nonlinear  spring  force  between  m*  and 
ff»i :  kiX  .  +  ptxi 

(v,)' 

scaled  velocity  at  1  =  t. 

X 

relative  displacement  between  mass  and 

A 

mean  value  of  the  modified  nonlinear 

foundation  for  one-degree-of -freedom 

spring  force  for  t  m<  t  <  , 

system 

h 

finite  time  increment;  At 

X' 

relative  displacement  between  mass  and 
foundation  at  t  =  tm  for  one-degree-of- 

k, 

spring  constant  between  m,  and  m,  , 

freedom  system 

m, 

ith  mass 

X, 

relative  displacement  between  mi  and 

W » -  1 

P 

«V  1  —  a* 

(*.)» 

relative  displacement  between  m.  and 
m,  i  at  /  *  tm 

r 

Vi  -  a1 

y 

absolute  displacement  of  mass  in  one- 

s. 

F'.,  -F' 

degree-of-freedom  system 

Si  ■ 

S'  —  S'  i  =  h  '  •  i  —  2f '  +  A '  i 

y- 

absolute  displacement  of  mass  at  t  —  tm 
for  one-degree-of -freedom  system 

s. 

Z ' ,  i  —  Z' 

y< 

absolute  displacement  of  m. 

(yi). 

altsoluie  displacement  of  mt  at  t  —  t. 

si  . 

S»  ~  S»  i  =  Z',i  2Z '  +  Zm  i 

z 

absolute  foundation  displacement 

i 

inde|K'ndent  time  varialile 

at 

r./2m,w. 

u 

scaled  velocity  i/w  for  une-degree-ol- 

coefficient  of  nonlinear  term  of  spring 

I  reedom  system 

y? 

ftilm, 

U, 

sc  aled  veloc  ity  at  1  =  1'  for  one-degree- 

ol-freedom  system 

ft. 

della  term  lor  the  ith  expiation  of  motion 

u, 

scaled  velocity  Xi/«i 

» 

w  k 

(U,)' 

scaled  veloc  ity  at  1  =  1, 

S 

k,lmi 
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A  Numerical  Method  for  the 
Transient  Response  of  Nonlinear  Systems 

C.  J.  O'Hara  and  P.  F.  Cunnid 


Structures  Branch 
Mechanics  Division 

Numriual  initiation  r<|ii.ilium  ait*  d«*nv«*<l  tot  drlri  mining  the  it-sjioiisc  of  iioiiIiikmi  systems 
siibjcdnl  lo  liaiisicnl  loads.  I  hr  minirntal  method  tonsisls  ot  .i)>|>ioxinidting  ihr  nonlinrai  vaiiahlrs 
and  ihr  lotting  (tint lions  in  dir  diflrirnital  equations  ovri  a  shoil  inlet \al  ul  uinr  by  dim  inr.in 
value,  by  a  straight  line,  ot  by  a  patalxtla.  I  bis  allows  lot  Dtihatnrl  integral  t\|>e  solutions  lot  ihr  non- 
linrat  trims.  A  step  by  step  solution  lollows  uhith  list's  an  itnalioii  method  timing  rath  iiuinnrnt  ol 
ihr  solulion.  I  hr  stillitirnt  ttmtliiitm  loi  ihr  tomrtgrntr  ol  the  itnalioii  method  is  prrsrntrd  bit  t hr 
tasr  ot  N  nutiiriit al  equations.  A  staling  law  is  pirsmlrd  which  rliininatrs  linrai  damping  tioin  thr 
rtpiation  ol  motion  by  a  pirst  i ilx*tl  itansloi  ination.  Kxamplr  piotilrms  ot  a  oticalegicc-nl-liecdom 
svsirni  anti  a  iwn-degiec-oMiecdoin  systrin  air  solved  by  thr  iiiiinrnt.il  mlrgi alion  equations  anti  the 
solutions  air  iniiipaird  vs  it  h  irsjxuisr  tutvrs  obtained  lioin  analtig  toiiiputns  at  \KI.. 


INTRODUCTION 

This  report  deals  primarily  with  approximate 
numerical  solutions  of  a  single  or  a  set  of  non- 
aiitonomous  second-order  ordinaiy  nonlinear  dif¬ 
ferential  equations.  While  the  c  lass  of  problems 
under  consideration  lie  in  the  field  of  structural 
dynamics,  the  promised  solutions  are  applicable 
to  many  other  physics  and  engineering  fields. 
The  mathematical  tools  of  ordinary  nonlinear  dif¬ 
ferential  equations  are  particularly  useful  when 
dealing  with  autonomous  solutions  or  with  ap¬ 
proximate  steady-stale*  solutions  of  these  prob¬ 
lems.  However,  they  require  considerable  ingenu¬ 
ity  and  insight  to  apply  and  are  not  suited  for 
the  study  of  transient  behavior.  An  attempt  is 
made  in  this  re|iort  to  present  an  easily  under¬ 
stood,  yet  |Miwerful  and  precise,  technique  which 
will  allow  most  engineers  to  cu|x*  with  the  tran¬ 
sient  res|ionse  of  nonlinear  systems. 

Two  previous  NRI.  Re|  torts  (1,2)  have  dealt 
with  this  problem,  and  this  rejiort  pursues  the 
same  general  approac  It.  Most  approximate  nu¬ 
meric  al  techniques  tail  lo  attack  directly  the  non¬ 
linear  dilferential  equations  in  their  solutions. 
Rather,  they  introduce  Mac  lauriu  or  Taylor  series 
exclusions  of  the  functions  as  in  the  method  of 
Picard  (.1).  The  othci  general  approach  is  lo  re¬ 
place  dillerential  equations  hv  equations  of  finite 

NRI.  Pudilrm  RO.V24R,  hnfi t  WW  (Ml  I  his  is  <«ti  fiitrtiin  irpmi 
im  nnr  ol  thr  |itnMrni.  tv  (ontmiimg 

M^ixivii|H  vilmiitlnl  Itnrmlirt  II,  I’lbJ 


differences  and  to  use  these  ecpiations  as  an  ap¬ 
proximation  to  the  differential  equations.  These 
are  good  general  pur|M>se  techniques.  However, 
they  lend  to  Ik-  routine  techniques  which  remove 
the  analyst  Irout  a  clear  understanding  of  the 
manner  in  which  the  differential  equations  were 
solved. 

The  numerical  method  presented  uses  only 
those  mathematical  tools  which  are  lamiliar  to 
most  engineering  graduates  and  are  applied  di¬ 
rectly  to  the  class  of  dilferential  equations  under 
study.  It  should  not  lx*  construed  that  this  is  a 
crude  technique  and  that  the  solutions  will  lx* 
greatly  in  error  or  will  have  inherent  instabilities 
of  large  magnitude.  The  examples  in  this  re|x>rt 
show  the  op|M>sile  to  be  true. 

To  those  readers  who  are  already  familiar  with 
Refs.  I  and  2.  this  rc|tori  is  a  direct  application 
of  the  principles  explained  therein.  For  those* 
|x*rsons  who  have  not  read  them,  however,  it  is 
noted  that  this  rc|xtrt  is  completely  self-contained 
and  these  references  are  not  required  reading. 

BACKGROUND  THEORY 
The  Linear  Problem 

It  will  lx*  Ix-nefic  ial  to  review  a  numerical  inte¬ 
gration  melhcxl  (I)  which  is  used  to  solve  linear 
single-dcgrce-of-freedom  problems  before  pro¬ 
ceeding  tc*  the*  nonlinear  ones.  Consider  the  un- 
dam|M*d  lineal  oscillator  shown  in  Fig.  I  subjec  t 
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t-itf.  1  -  Ijncdi  miilljloi 

T 

r  (♦) 


lo  an  applied  Imre  F(t).  let  y  be  the  displace¬ 
ment  ol  the  oscillator  so  that  the  diHeienti.il 
equation  of  motion  is 

my+ky  =  F(t)  (I) 

where  the  dots  denote  differentiation  with  respect 
to  time.  II  oi 1  =  k/m,  tq.  (I)  Itetomes 

•  .  *  *’(<) 

y  4-  atty  =  — — .  (2) 

rfl 

The  general  solution  of  this  equation  is 

y  =  yr  +  »  (S) 

where  yr  is  the  complementary  .solution  and  yP 
is  the  particular  solution,  This  pro|>erty  of  linear 
equations  will  lie  shown  to  have  some  value  in  the 
approximate  solution  of  nonlinear  equations. 

For  the  case  under  study  the  complementary 
solution  is  well  known  and  the  particular  solu- 
is  a  Duhamel  integral  (4).  If  v  =  y/ot,  the  general 
solution  of  Eq.  (2)  is 

y  —  y0  cos  otl  +  Vu  sin  oil 


into  equal  segments  of  time/  and  represented  in 
some  approximate  manner  for  each  increment,  a 
step  by  step  approximate  solution  follows.  It  is 
noted  that  F.q.  (4)  is  true  for  all  times  during  the 
tes|M>nse  of  the  oscillator.  For  example,  suppose 
y  =  yi  and  v  =  v<  at  /  =  1 1.  Now  the  time  can  be 
redefined  arbitrarily  to  start  at  zero  for  the  next 
inc  rement  with  y >  and  Vi  being  the  initial  condi¬ 
tions.  I  he  Duhamel  integrals  are  solved  for  this 
next  step.  yt  and  vt  are  found,  and  a  repetition 
of  the  process  defines  the  next  pair  of  points. 
The  process  is  self-perpetuating. 

The  problem  in  this  direct  attack  upon  the  dif¬ 
ferential  equation  of  motion  has  resolved  itself 
into  the  solution  of  these  integrals  for  a  short 
time  increment.  Since  the  forcing  function  may 
he  known  only  as  a  graphical  function,  as  a  dis¬ 
continuous  function,  or  as  a  complicated  analytic 
function,  some  methods  of  describing  it  over  the 
immediate  range  of  integration  is  now  discussed. 

Approximate  Methods  of 
Representing  Functions 

Three  methods  are  presented  for  the  approxi¬ 
mate  representation  of  a  function  over  finite  in¬ 
crements  of  lime  A l  —  h.  Suppose  a  portion  of  an 
arbitrary  function  F(t)  is  divided  into  equal  seg¬ 
ments  of  time.*  Figure  2  shows  the  rectangular 
step  representation  consisting  of  horizontal  lines 
drawn  through  the  mean  value  of  the  function 
over  each  increment  h.  Appendix  A  reviews  two 
common  procedures  for  obtaining  graphically  the 
mean  value  of  a  function  for  a  given  increment. 
The  equation  of  the  function  during  each  incre¬ 
ment  is 

F{t)  =  F*  =  constant,  f»  <  I  <  f«.i  (5) 


+  — f  F(T)  sin  o)(t  —  T)  dT  (4a) 
mot  J, 

and  its  scaled  derivative  v  is 

v  -  -yo  sin  cot  +  t>»  cos  cot 

+  — f  F(T)  cos  oi(t  —  T)  dT  (4b) 
mot  Jt 

where  y,  and  v«  are  the  initial  values  at  t  =  t,  =  0. 

Usually  the  integrals  of  F.q.  (4)  cannot  be  eval¬ 
uated  for  an  arbitrary  curve  of  F.  If  F  is  divided 


where  Fm  is  the  mean  value  of  F  during  the  in¬ 
crement. 

The  second  method  represents  the  curve  by  a 
straight  line  through  the  end  (minis  across  each 
increment  as  shown  in  Fig.  .1.  The  equation  of  the 
(line  lion  from  t »  to  tn ,  i  is 

F(t)  =  F,  +  SK^  (6) 

where  .S’«  =  Fm* i  ~  F»,  and  time  begins  at  f„. 

‘Dividing  F  min  equal  tegmenta  of  lime  n  an  unnrtnvirv  but  ttm- 
venient  rrunt  ik hi  mint  M  make*  t  alt  ulatinm  easier  and  rom|Hiier  pnr- 
gramming  less  <  iimhertnme. 
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Fig  2  -  Krprrwnuiiiin  ol  .i  I uiu (ion  by  iccUiiKiilar  stt*|>s 


The  third  method  of  approximating  a  I  tint' lion 
across  an  increment  is  to  pass  a  parabola  through 
three  successive  |>oinls.  For  example.  Fig.  3  shows 
three  successive  |>oints  on  the  curve,  namely, 
F»  i,  A’«,  and  A'».  ■.  The  equation  ol  the  parabola 
from  t,  to  f.,i  passing  through  these  three  |M>inls 
on  the  curve  is 


F(t)  =  A, 


-(A 


') *A.  ,  +  [(h')*-h*]F'  +  h*FK.>) 

- **'(*  +  * t - r 


—  (A  -f  h')F ,  +  hF».  i 
AA'(A  +  A') 


<10> 


or 


r<o-r.*s.i  +  3(£-i)  a.„ 

r„)-f.  +  s.i  +  5U(  11 -i)  ,7b, 


Solution  Equationi  for  the 
Undamped  Linear  Oscillator 


If  the  parabolic  representation  of  a  function 
given  by  Eq.  (7b)  is  substituted  for  F  in  Eqs.  (4), 
and  the  integrations  performed  and  evaluated  at 
t  =  A,  there  results 


where 

SI  =  S, . ,  -  ,S.  =  Fm .  c  -  2A’, ,  i  +  A’.  (8) 

.S*  ,  =  Sm  -  Sm  ,  =  A', .  i  -  2Fm  +  A’.  ,  (9) 


i  =  y n  cos  9  +  vm  sin  9 

.S» /,  sin  9\ 

+  T  (i  ~  co*  +T  (i — T) 

SJ-,  [sin*  _  2(1  -  cos  g)l 

*2A  I  9  9*  J 


(11a) 


and  lime  again  begins  at  t„.  Note  that  .S’  is  not  the 
square  of  S.  Equation  (7a)  is  used  to  represent  the 
curve  during  the  lirst  increment  of  time,  that  is, 
lor  n  ir  0.  Equation  (7b)  is  the  expression  used  for 
the  remaining  segments  of  the  curve.  Equation 
(7a)  is  derived  s|>ccially  to  avoid  the  nonexistent 
let m  in  .S'*,  (that  is,  F  ,)  if  Eq.  (7b)  were  used  for 
the  Inst  increment. 

for  the  c  ase  where  two  sun  essive  inc  rements  are 
not  equal,  say  A  between  A',  t  and  A',  and  A' 
lirtween  A'„  and  A',,i.  Eq.  (7b)  is  adjusted  to  read 


fm  i  =  —y*  sin  1  +  r,  cos  9 

+  £,ing  +  s.L>-j™  *1 

k  k  9 

^(L m« 

where  9  =  wA.  It  is  noted  that  for  constant  incre¬ 
ments  A,  the  trigonometric  coefficients  are  calcu¬ 
lated  just  once  for  the  entire  solution. 
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Equations  (11)  are  used  for  the  straight  line 
representation  of  F  by  setting  S*  equal  to  zero, 
and  for  the  rectangular  step  representation  by 
setting  both  S  and  S*  equal  to  2ero.  The  equa¬ 
tions  for  the  latter  case  are 

y**i  =  y»  cos  0  +  vm  sin  0 
V 

+  ^0  “«>*«)  (12a) 


t»,n  =  —y,  sin  0  +  vm  cos  0 


sin  0. 


(12b) 


Recall  that  is  the  mean  value  of  the  function 
during  each  increment.  For  functions  such  as 
F,  which  are  explicit  functions  of  time,  it  is  recom¬ 
mended  to  find  the  average  value  of  F  during  the 
time  increment  instead  of  employing  the  graphical 
techniques  outlined  in  Appendix  A.  The  average 
value  of  the  function  represented  by  a  straight 
line  is 


F.  «  -  -"  ‘-1  (13) 

2 

while  the  average  value  for  a  parabolic  represen¬ 
tation  is 


F.  “  yjj  (~F»  t  +  *F,  +  5F\,,i).  (14) 

hither  F.q.  (13)  or  (14)  is  then  substituted  for  Fm 
in  F.qs.  (12)  and  the  numerical  integral  equations 
solved  in  a  step  by  step  fashion. 

Solution  Equations  for  the 
Viscously  Damped  Linear  Oscillator 

The  method  of  solution  for  the  undauqted  case 
is  also  applicable  to  the  viscously  damped  sys¬ 
tem.  The  equation  of  motion  of  an  oscillator  with 
linear  damping  is 

my  +  cy  +  iry  =  F’(<) .  (15) 

If  a  =  r/2mco,  F.q.  (15)  becomes 

y  +  2  nmy  +  cu*y  =  tAil  (16) 


The  form  of  the  general  solution  for  Eq.  (16) 
is  given  by  Eq.  (3)  which,  for  this  case,  consists 
of  a  combination  of  exponential,  trigonometric, 
and  hyperbolic  functions,  and  the  solutions  may 
be  found  in  Appendix  B.  The  forcing  function 
is  approximated  as  previously  discussed  and  the 
solution  carried  on  precisely  as  in  the  undamped 
case.  For  example,  with  a  =  1  and  the  average  or 
mean  value  of  the  force  used,  the  equations  are 

y.4 1  =  y»(l  +  0)e~*  +  v»0e 

+  ^  [1  -  (1  +•)«-•] 

V 

Vm  *  i  =  -y.  0t~*  +  v»(  1  —  0)e~*  +  0e~*. 


Comments 

The  solution  of  the  linear  problem  leads  to 
several  interesting  observations. 

1.  A  direct  attack  on  the  differential  equation 
is  made. 

2.  The  coefficients  of  the  variable  terms  in  the 
numerical  equations  become  constants  through¬ 
out  the  solution  when  equal  time  intervals  are 
employed. 

3.  The  solution  makes  use  of  the  natural  expan¬ 
sion  functions  for  the  differential  equations.  That 
is,  they  are  of  the  form  of  trigonometric,  hyper¬ 
bolic,  and  exponential  f  unctions  as  they  would  be 
in  an  analytical  case. 


NONUNEAR  PROBLEM 
Numerical  Solution 

Suppose  an  ordinary  second-order  differential 
equation  is  reducible  to  the  form 

y  +  H{y,  y,  »)  =  0  (17) 

where  H[y,  y,  l)  may  be  a  very  complicated  func¬ 
tion.  Since  the  function  H  contains  a  forcing  func¬ 
tion  F,  Eq.  (17)  bec  omes 


y 


where  -»  <  a  <  «,  and  w  >  0. 


(18) 
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r» 


where 


H  =  C 


F 

m ' 


I  .el  K(|.  (18)  Ih'  rewritten 


numerical  integration  equations  for  y  and  v  are 
obtained.  Consider  a  simple  osc  illator  with  a  cubic 
hardening  spring.  The  equation  of  motion  is 

my  +  ky  +  fly3  =  F(t) 


or 


y  +  2aou>  +  coJy  =  — 
m 


y  +  oj*y  = 


Fit) 

m 


~[G  —  2auy  —  u*y] 

and  let 

wJ8  =  6'  —  2a(uy  —  tu-y 


where 

6=f  y*- 


so  t  hat 


Use  Kq.  (7b)  for  the  cubic  term,  obtaining 


F 

y  +  2auy  +  ca’y  = - 01*8 .  (1 9a) 

m 

If  8  were  zero  for  all  time,  the  equation  would 
be  linear,  and  the  solution  has  already  been  pre¬ 
sented.  For  convenience  let  the  damping  term 
a  Ih *  zero  in  F.q.  ( 19a).  Then, 

F 

y  +  cu’y  = - ci»*8 .  ( 1 9b) 

m 

The  solution  of  this  equation  is  similar  to  Kqs. 
(4);  that  is, 

y  —  y«  cos  8  +  vo  sin  0  +  —  f  F(T) 

mu  l 


a-f  [y.+  (y».i  ~  y»)  j 


(21) 


F.xpand  Fq.  (21),  substitute  it  into  Kqs.  (20)  for 
8,  and  integrate  each  term  over  the  increment. 
If  the  forcing  function  were  also  approximated 
by  the  parabolic  representation,  the  resulting 
numerical  integration  equations  are  of  the  form 
of  Kqs.  (II)  with  the  additional  terms  ym.t,  y«. 
and  ya .  1 ;  that  is, 


sin  u(l  —  T)dT 


-r 


8(D  sin  w(t  -  T)dT 
(20a) 


y»o  1,  Fit,  F H4i,  y.-i,  y»,  y»»i,  ti»,A) 

(22a) 


v  =  —ya  sin  8  +  va  cos  8  4- 


-  r  * 

mu  l 


F(T) 


e»«i  gi  ( F  „  - 1 ,  F a,  F ■ ,  1 ,  y«-i,  y»,  y»*  1 ,  vn,h ) . 

(22b) 


cos  u( t  —  T)dT  —  c a  f  8(T)  cos  u(t  —  T)dT. 

%  (20b) 

The  solution  of  the  linear  problem  for  a  given 
arbitrary  curve  of  F  requires  that  F  be  partitioned 
into  finite  increments  of  time  and  be*  approxi¬ 
mated  over  each  increment  by  one  of  several 
representations.  The  same  approach,  using  the 
same  increment  h ,  is  prcqioscd  to  handle  the  inte¬ 
grals  in  Kqs.  (20)  containing  the  nonlinear  terms 
in  8. 

For  example,  each  term  in  8  might  he  approxi¬ 
mated  by  the  fiarabolk  representation  so  that 


Kvcrylhing  is  known  on  the  right-hand  side  of 
Kqs.  (22)  except  y„,  1.  As  a  first  trial  this  value  may 
Ih-  assumed' equal  to  y„.  Substitute  this  into  the 
right-hand  side  of  Kq.  (22a)  to  find  y»»  t.  Use  this 
value  for  y„ ,  ■  in  the  right-hand  side  of  Kq.  (22a) 
to  find  a  second  value  of  y„,  (.  Repeat  this  itera¬ 
tion  process  until  the  succeeding  values  of  y„,i 
converge.  Use  the  final  value  of  y„,i  in  Kq.  (22b) 
to  find  Cut  1. 

T  his  method  of  solution  could  also  be  used  if 
instead  of  the  parabolic  representation  the  straight 
line  representation  approximated  the  nonlinear 
term.  However,  in  either  c  ase  a  great  deal  of  time 
and  effort  is  required  if  the  analyst  uses  a  desk 
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calculator.  Of  course,  the  numerical  equations 
could  be  programmed  for  an  electronic  computer. 
Everything  which  follows  from  this  point  is  di¬ 
rected  toward  desk  calculator  computations. 

The  recommended  method  of  solution  with  a 
desk  calculator  is  to  use  the  mean  value  or  average 
value  of  the  variables  in  the  &  term  during  a  given 
finite  increment.  If  this  assumption  is  also  ex¬ 
tended  to  the  forcing  func  tion,  Eqs.  (20)  are  inte¬ 
grated  to  give 

y* ,  i  =  y»  cos  8  +  vm  sin  8 

+  4-*  ( 1  -  cos  8)  -  S.<  1  -  cos  8)  (23a) 

A 

v, .  i  =  — y„  sin  8  +  v,  cos  8 

+  sin  8  —  5*  sin  8.  (23b) 

k 

For  the  oscillator  with  the  cubic  hardening  spring, 
the  straight  line  averaging  method  yields 

6.  =  ^  (y.  +  y...)3  (24) 

while  the  parabolic  averaging  method  yields 

=Jj28  * ~y "  1  +  8y»  +  5y»*i)s-  (25) 

Once  again  y« ,  ■  is  the  unknown  cm  the  right- 
hand  side  of  Eqs.  (23)  and  the  iteration  process 
is  used  for  each  step  of  the  solution  as  previously 
described. 

The  First  Increment 

Special  treatment  is  necessary  for  the  applica¬ 
tion  of  the  numerical  integration  equations  at  the 
start  of  the  solution.  At  the  beginning  of  a  problem 
the  initial  conditions  are  always  known  so  that 
y,  and  t/«  are  established.  Having  selected  a  time 
inc cement  h,  the  first  step  of  the  solution  depends 
upon  the  type  of  averaging  method  to  be  used  for 
the  nonlinear  terms.  In  the  case  of  the  straight 
line  averaging  method,  the  first  trial  value  for 
y(  might  be  assumed;  yi  might  be  set  equal  to 
y«,  or  a  Mac  laurin  series  might  be  used  to  approxi¬ 
mate  yi.  In  any  event,  the  first  trial  value  of  y,  is 
substituted  into  the  right-hand  side  of  the  nu¬ 
merical  integration  equation  to  find  a  new  value 
of  yi.  The  iteration  procedure  follows  as  previ¬ 


ously  mentioned.  Of  course,  if  5  contains  the 
scaled  velexity,  the  same  approach  is  used  to  find 

Vf 

In  the  case  of  the  parabolic  averaging  meth¬ 
od,  Eq.  (25)  shows  the  rear  point  y»- 1,  which  must 
be  known  before  iterating  to  find  y**i.  At  the 
start  of  a  problem  where  n  =  0,  y  - 1  does  not  exist. 
It  is  suggested  that  the  straight  line  average  be 
used  for  the  first  increment  to  establish  y,.  If 
greater  accuracy  is  desired,  use  the  linear  averag¬ 
ing  method  for  half  an  increment,  that  is,  for 
A/2.  Having  established  y m  in  this  manner,  the 
parabolic  averaging  method  is  now  used  for 
another  half  increment  to  find  yi.  The  full  in¬ 
crement  might  then  be  used  from  this  point 
throughout  the  remainder  of  the  solution. 

Graphical  Form  of  Nonlinear  Components 

Quite  frequently  the  nonlinear  characteristic 
of  a  material  in  a  system  is  determined  from  lab¬ 
oratory  experiments  and  is  plotted  as  force  ver¬ 
sus  displacement  or  velexity.  It  is  sometimes  pos¬ 
sible  to  find  an  analytical  expression  for  such  a 
curve.  In  the  event  this  is  not  readily  attainable, 
a  graphical  technique  for  finding  the  mean  value 
of  a  function  over  an  increment  is  used  (see  Ap- 
pendix  A). 

For  purposes  of  illustration  consider  Fig.  4a, 
which  shows  the  spring  force  for  positive  dis¬ 
placements  only.  A  line  tangent  to  the  curve  at 
the  origin  is  drawn  and  is  labeled  the  A-line.  Fig¬ 
ure  4b  represents  the  spring  force  minus  the 
A-line  as  a  func  tion  of  displacement.  This  is  the 
curve  from  which  the  mean  value  of /is  deter¬ 
mined.  Suppose  an  oscillator  contains  such  a 
spring.  Equations  (23)  become 

y»« i  =  y»  cc>s  8  +  vn  sin  0 

+  ^  ( 1  —  cos  0)  —  ^  ( 1  —  cos  8)  (26a) 

i  =  — y»  sin  8  +  v.  cos  8 

+  y  sin  8  -  ^  sin  8  (26b) 

where  F»  is  the  mean  value  of  the  forcing  function 
and  f,  is  the  mean  value  of  the  curve  shown  in  Fig. 
4b.  At  step  n,  y»  and  v *  are  known,  F„  is  known 
from  the  input  curve,/*  depends  upon  y»  and  the 


U.S.  NAVAL  RESEARCH  LABORATORY 


7 


(a)  Force  displacement  curve 


Fix  4  -  Spring  lore  col  a  nonlinear  spring.  Ihc- time  is  tangent 
to  the-  <  urve  at  the  origin. 


unknown  y«t|.  As  a  first  trial  for  finding  y„t|, 
find  /  for  y„  and  use  this  for  /„.  Substitute  into 
Eq.  (26a)  to  find  a  first  trial  of  y„«|.  Now  find 
fn  between  y,  and  y» , i  from  Fig.  4b  and  substitute 
this  into  Eq.  (26a)  to  obtain  a  new  value  of  yHll. 
Kc|x*al  the  process  until  succeeding  values  of 
y*  *  i  converge. 

For  certain  curve  sha|>es  it  might  lie  advan¬ 
tageous  to  draw  two  or  more  A-lines  for  certain 
|Kirtions  ol  a  given  curve.  These  A-lines  would 
provide  a  solution  which  follows  more  closely  a 


piecewise  linear  solution  of  the  system  by  reduc¬ 
ing  the  magnitude  of  the  adjusted  forces /.  This 
means  that  a  corresponding  number  of  differen¬ 
tial  equations  must  lie  written  for  each  region  of 
the  curve  where  a  A-line  is  drawn.  Proper  initial 
conditions  and  cu's  must  lie  determined  for  each 
numerical  integration  equation.  An  example 
might  be  a  material  whose  force-displacement 
curve  follows  closely  an  ideal  elastic-plastic  re¬ 
lationship. 

Forcing  Functions 

Foundation  motion  of  structures  is  an  important 
type  of  forcing  function  in  the  field  of  structural 
dynamics.  Such  motion  may  be  described  as  foun¬ 
dation  acceleration,  velocity,  or  displacement.  In 
the  case  of  foundation  acceleration  the  differen¬ 
tial  equation  of  motion  for  a  nonlinear  oscillator 
with  a  cubic  hardening  spring  is 

*  +  *>**  =  -  i'-y**3  (27) 

where  x  is  the  relative  displacement  between  the 
mass  and  the  foundation  and  i‘  is  the  foundation 
acceleration.  This  equation  is  similar  to  Eq.  (2), 
with  x  replacing  y  and  ~i’  replacing  F/m.  The 
parabolic  averaging  method  is  recommended  for 
systems  with  a  known  curve  for  a  This  average 
lor  i  is  the  same  as  Eq.  (14)  provided  the  F  terms 
are  replaced  by  i‘  terms.  Equations  (II)  may  be 
also  used  provided  the  following  changes  are 
made: 


•S*.  SI., 

a  w*  *  y*  •r"’ v*  ~  **"• 

When  the  foundation  velocity  is  the  prescribed 
input,  an  interesting  relationship  is  found  for 
the  parabolic  average  of  z  to  be  used  in  the  nu¬ 
merical  solution  of  Eq.  (27).  Consider  the  para¬ 
bolic  representation  for  foundation  velocity  and 
differentiate  to  find  foundation  acceleration: 
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The  average  acceleration  during  the  increment  is 


This  is  the  approximate  representation  lor  the 
foundation  acceleration  according  to  the  theorem 
of  the  mean  of  differential  calculus. 

II  the  Inundation  displacement  is  the  prescribed 
input,  the  differential  ccjuation  of  motion  for  the 
nonlinear  oscillator  should  lie  of  the  form 

yt  +  tidy,  =  adz  -  y*(y,  -  r)s  (28) 

where  y%  is  the  absolute  displacement  of  the  mass. 
I'he  paraltolic  average  representation  for  z  is 
lotincl  for  eac  h  increment  from  the  given  c  urve. 

Step  Change*  in  Input 

l  he  numerical  integration  equations  are  par¬ 
ticularly  applicable  to  problems  in  which  step 
changes  occur  in  the  input  of  the  system.  The 
only  restriction  in  the  solution  is  to  require  that 
.1  given  increment  terminate  where  the  step  change 
occurs,  let  the  curve  shown  in  Fig.  5  Ik*  a  foun¬ 
dation  velocity  of  a  nonlinear  osc  illator,  and  sup¬ 
pose  the  increment  h  is  selected  for  the  solution. 
After  the  fifth  step  of  the  solution,  a  new  incre¬ 
ment  h'  is  used  for  the  sixth  inclement  in  order 
to  arrive  at  the  time  when  the  step  c  hange  in  the 
foundation  velocity  occurs.  The  si/e  increment 
h  is  then  used  throughout  the  remainder  of  the 
solution. 


Foi  the  sixth  increment  of  the  numerical  inle- 
gial  solution,  a  new  increment  O'  =  c oh'  must 
if  place  0  =  i ah  in  the  solution  equations.  Equa¬ 
tion  (10)  must  Ik*  used  it  the  input  is  being  ap- 
pioxitnated  by  the  paralmlic  representation.  In 
so  doing,  the  solution  is  found  at  step  6  for 
and  Uh,  where  xK  is  the  relative  displacement  be¬ 
tween  the  mass  of  the  oscillator  and  the  founda¬ 
tion  and  <c«  is  the  corresponding  scaled  relative 
velcKity.  At  this  |M>int  the  step  change  in  foun¬ 
dation  veliHity  (scaled)  must  Ik*  added  to  u«  to 
give  u».  With  xH  and  u»  as  the  initial  conditions 
for  step  7  of  the  solution,  and  using  0  in  the  solu¬ 
tion  equations,  the  input  curve  should  be  repre¬ 
sented  by  the  straight  line  method  between  points 
6  and  7  instead  of  the  parabolic  methcKi  due  to  the 
discontinuity  at  point  H.  I  bis  completes  the  solu¬ 
tion  to  step  7.  The  parabolic  representation  may 
Ik*  used  from  this  |K>int  to  the  end  of  the  input 
c  urve*. 

CONVERGENCE  OF  THE 
ITERATION  METHOD 

One  may  wonder  if  an  assumed  variable  could 
not  produce  a  diverging  solution  when  substi¬ 
tuted  into  the  5  term  on  the  right-hand  side  of 
the  numerical  equations  for  a  given  increment. 
Scarborough  (3b)  has  discussed  the  sufficient 
conditions  for  convergence  of  one  and  two  nu¬ 
merical  equations  when  the  iteration  prexess  is 
applied.  An  extension  of  the  pre wed ure  follows 
for  a  set  of  A  numerical  ecpiations. 

Consider  the  set  of  equations 

c,  =  *2,  *.%],  i  =  l,2 . A  (29) 

where  B ,  is  a  set  of  known  functions  in  terms  of 
the  c/s.  The  set  of  ecpiations  is  satisfied  by  the 
exact  values  of  the  set  of  roots  cj.  As  a  first  ap¬ 
proximation  to  find  a  set  of  roots,  try  c1*'  ■  Hence, 
Eq.  (29)  gives 

«/"  «*,W\  «r . «H*  (SO) 

Substrata  Eq.  (30)  from  (29)  and  apply  the  theorem 
of  the  mean  value  for  a  function  of  A  variables. 
This  gives 

b,j 

j-i 


tig.  r>  Slc-|>  (hatiKT  in  inuiicbimn  velocity 


(.31) 
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+  »(€,  ~  Q . t(s01  +  cfrU.V  ~  cj°*)  ] 

dtj 

0  *£  <f>  S  1. 

Add  each  equation  of  the  set  of  equations  as  ex¬ 
pressed  in  Eq.  (31)  and  consider  only  the  absolute 
values.  Thus, 

t  £  t  Uj-«i0,ll*J-  <»*> 

i  i  <- 1  i- 1 

Let  the  maximum  value  of  the  terms  {|6i,i|  +  ... 

+  |«X.,U . {|Bi..v|  +  ...  +  |0v.*|}  l>e  a  proper 

traction  <ft  tor  all  points  in  the  region  <,). 

Then  Eq.  (32)  becomes 

2  <»s> 


This  relation  holds  for  the  first  approximation. 
For  succeeding  approximations,  similar  relations 
are  obtained.  That  is. 


2  k  -  *  2  k  -  «,"l 


as  desired  t>y  rejR-ating  the  iteration  pnxess  a 
sufficient  number  of  times.  This  means  that  the 
errors  | e,  —  e,<M|  can  be  made  as  small  as  desired. 
Therefore,  the  iteration  process  converges  when 
the  N  conditions 

IJT...I  +  |0*., I  +  •••  +  |0.v.,|  <  1 


|0|,.v|  +  |0*.«|  +  •••  +  |0.V,A'|  <  1  (36) 

are  satisfied  for  all  |x>ints  in  the  neighhorhcxxl  of 

«r 

Consider  Eqs.  (23),  which  pertain  to  an  incre¬ 
ment  of  time  h  for  a  single-degree-of-freedom 
system.  All  terms  on  the  right-hand  side  of  the 
equations  are  constants  except  5*.  Suppose 
is  a  function  of  displacement  y  and  scaled  ve¬ 
locity  u.  Equations  (23)  are  of  the  form  of  Eq. 
(29)  where 

~  y*  *  i  *  *2  ~  ♦  i 

01  =  yn  cos  0  +  Vn  sin  6 


+  ( 1  -  ‘os  9)  —  &,(  1  —  cos  9) 

Bi  =  —  y»  sin  9  +  v,  cos  0 
f 

+  -j*  sin  9  —  &>  sin  9. 


2  « #2  i«* - «/*  "I-  <*»> 


According  to  Eq.  (36)  the  iteration  pnxess  on  the 
fin's  converge  in  the  nth  step  of  the  solution  pro¬ 
vided  that 


Multiply  together  all  these  inequalities  as  ex¬ 
pressed  in  Eqs.  (33)  and  (34)  and  divide  through 
by  the  common  factors 

2  k  -  2 1‘‘  ~  «‘,,l . 2  k-«i*'nl 


( 1  —  cos  0) 


+  sin  0  ~~  <  I 

fiy..  i 


( 1  —  cos  0)  +  sin  0  -  <  1 

OVul  Ot'MI 


in  the  neighborhexxi  of  [ y',0,1 , ,  r,„04,l  ] . 


so  that. 


2  k-«H  ^  2  k-«Swl-  (3S) 

C-l  t-t 

Since  f  is  a  proper  function,  the  right-hand 
ntemlxT  of  this  inequality  may  lx*  taken  as  small 


A  SCALING  LAW 

Outsider  the  equation  of  motion  for  a  single- 
degree-of-freedom  system: 

y  +  2  auty  +  miy  =  — m*8  (37) 
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(a)  System 


(b)  Foundation  input 


U)  *  \s  I  obtamrd  from  the  numerital  method  and 

(mm  an  analog  computet 


(■»K  •*  *  Input  and  numerical  solution  o(  a  nonlinear  singlc-dcgrer-oMrerdom 
s* st fin  with  vis< ous  damping 


where  8  =  8(y, y,  < ).  Substitute  the  transformation 
y  =  Ae"“'  (38) 

into  K,l  (37).  'I'his  gives 

A  +  p*A  =  -w’ip"*'  (39) 

where  p*  =  w*  ( I  —  a* )  and  A  =  A  ( A ,  A ,  t ) . 

Equation  (39)  ran  Ik-  solved  graphically  (2)  or 
numerically  and,  using  the  transformation  of 
Ec|.  (38),  the  response  of  the  original  system  is 
found.  Fot  the  rase  where  the  mean  or  average 
value  of  (lie  terms  in  the  A  expression  are  used, 
the  numerical  integration  equations  lor  Kq.  (39) 
are 


Am,  i  =  A ,,  ros  pk  H — 5  sin  ph 
P 


—  A„  +  (1  —  ros  ph)  (40a) 

^  — -  =  —A»  sin  ph  +  —  ros  ph 
P  P 

—  A„  sin  ph.  (40b) 


EXAMPLES 
Example  I 

Figure  (ia  shows  the  single-degree-of -freedom 
system  subjected  to  a  transient  foundation  velocity 
shown  in  Fig.  Ob.  The  equation  of  motion  is 
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x  +  2ttw«  +  oj2x  =  — <11*6,  a  <  1 

where  8  =  fix3lk  +  i'/o)1.  Since  a  desk  calculator 
is  used  to  find  the  res|>onse  x  of  the  systent,  the 
mean  value  of  the  terms  in  the  6  expression  is 
used.  For  this  case  a  <  1.  so  that  the  equations 
from  Appendix  B  for  Case  II  apply.  By  replacing 
F/m  by  —ott8  in  F.q.  (Bl),  the  numerical  integra¬ 
tion  equations  are 

x,.i  =  x,  f  at  ^cos  ph  +  j  sin  ph^j 

e-a» 

+  u, - sin  ph 

r 

—  8m  [1  —  t  at  (cos  ph  +  j  sin  ph)]  (41a) 
e"a* 

Ur. i  =  -*« - sin  ph 

r 

(01  \ 

cos  ph  — -  sin  ph J  —  8m  — sin  ph. 

(41b) 

The  values  of  the  parameters  are 

m  =  0.10  Ih-secVin.  c  =  2.4  Ib-sec/in. 

k  =  1440  Ib/in.  0  =  720  lb/in.3 

A  time  increment  h  =  0,0028  sec  is  selec  ted  so 

that  8  =  360  wA/2ir  =  19.25  degrees.  Referring 
to  Fig.  6b  it  is  noted  that  twelve  steps  of  this  in¬ 
crement  arrive  at  f  =  0.0336  set  and  that  the  thir¬ 
teenth  step  requires  an  h'  =  0.0014  set  to  com¬ 
plete  the  input  of  the  system.  The  solution  of 
the  free  vibrations  from  this  point  is  based  u|>on 
h  =  0.0028  sec. 

Since  the  foundation  vekxity  is  the  known  in¬ 
put,  while  the  8  expression  calls  for  the  founda¬ 
tion  acceleration,  the  latter  is  approximated  over 
each  short  increment  by  the  parabolic  average 
tneth<»d.  That  is, 

••  in  *  I  —  *R 

*.  =  — Jj — . 

The  parabolic  averaging  method  is  used  for  the 
x3  term  in  the  8  expression.  Upon  substituting 
the  values  of  the  parameters  into  Fqs.  (41)  for 
h  =  0.0028  set ,  there  results 


(41a') 


Um .  i  =  -0.318880  Xn  +  0.881534  u„ 

-  0.318880  6,  (41b') 

where 

.  _  X, .  I  ~  Zn 

*  40.320 

+  0.000289  (— x, - 1  +  8  x,  +  5  x,.,)3.(41c') 


The  straight  line  averaging  method  is  used 
for  the  first  increment  of  the  solution  since  x  i 
does  not  exist.  Table  I  shows  the  arrangement 
for  solving  the  equations  and  presents  the  data 
for  n  =  0,  1,  and  2.  Although  all  numbers  are 
carried  to  six  decimal  places,  the  iteration  of  the 
variable  x  was  carried  to  four  decimal  places. 
This  generally  required  four  trials  in  each  step. 
The  results  of  the  numerical  method  are  plotted 
in  comparison  with  the  response  curve  obtained 
from  an  analog  computer  at  NRL  as  shown  in 
Fig.  6c. 


Example  2 

Figure  7a  shows  a  two-degree-of-freedoni  sys¬ 
tem  subjected  to  a  transient  foundation  accelera¬ 
tion  shown  in  Fig.  7b.  There  are  two  cubic  harden¬ 
ing  springs  in  the  system  with  the  following 
force-displacement  relationships: 

ft*  =  k\X\  +  0ixf,  /n  =  AjXj  +  0jxf. 

The  equations  of  motion  are 

it  +  2aio*iXi  +  w*xi  =  —  «i^8i  (42a) 

xi  +  2aj<mx2  +  u\xt  =  — »?8 1  (42b) 


where 


C|X|  _  w? 
miwf  oil 


x,.,  =  0.945311  *,  +  0.318880  u. 

-  0.054689  8» 
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lution  <>l  a  nonlinear  two-degree- 
of-freedom  system  with  viscous 
damping 
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The  numerical  integration  equations  for  Flqs. 
(42)  are  the  same  as  in  Example  1  with  the  addi¬ 
tional  requirement  that  the  proper  subscripts  In- 
used .  The  values  of  the  parameters  are 

nii  =  0.40  lh-sec*/in.  m-  =  0.10  lb-sec  Vin. 

k,  =  8000  Ib/in.  kt  =  2500  Ib/in. 

(i,  =  7000  lh/in.s  /3,  =  2000  Ib/in.® 

ri  =  12  ll>-sec/in.  r*  =  .'1  Ib-sec/in.® 

A  time  increment  h  =  0.002  sec  is  selec  ted  so 
that  0,  =  .460  cu,A/2ir  =  10.21  degrees  and  0,  = 
360  (u-hlin  =  18.12  degrees.  The-  paralmlic  av¬ 
eraging  methcKl  is  used  for  the  variables  in  the¬ 
ft  expressions  (including  the  input.  Fig.  7I>).  1'|m>ii 
substituting  the  values  of  the  parameters  into 
the  numerical  integration  equations,  there  results 

( x, ) » . ,  =  0.961066  <x,)„  +  0.270884  (u,). 

-  0.038934  <  ft,  )„ 

(u,)..,  -  -0.270884  (x,  >,  +  0.903603  (u,)„ 

-  0.270884  (ft,), 

where 

(  ft,  ).  =  0.875000  -  0.059293  u,  -  0.312500  x, 

-  0  250000  *  +  20 ^00 

and 

( x, ), . i  0.951391  (x.lff  )  0.301839  (u, )„ 

-  0.048609  (ft,), 

(«,)„.,  =  -0.301839  (x,)„  +  0.894121  («,), 

-  0.301839  (ft,), 

where 

(ft, »,  =  0.250000  v,  f  I  0.047434  «, 
-0.169706  u,  -  0.8(KMHM)  t,  -0.700000  r?. 


The  bars  altove  the  variables  in  the  6  expres¬ 
sions  i  epic-sent  paraltolic  average  values  during 
the  inclement  n  in  ft  +  I.  I  he  convergence  of 
die  variables  of  (lie  iteration  method  was  cairied 
to  four  decimal  places.  Ibis  generally  required 
five-  trials  in  each  step.  The  numerical  results 
are  plotted  against  the  response  c  urves  obtained 
horn  an  analog  computer  at  XRL  and  are  shown 
in  Figs.  7c  and  7d. 

CONCLUSIONS 

A  numerical  integration  method  has  been  pre¬ 
sented  whic  h  is  easily  understood  and  provides 
a  good  solution  for  the  transient  response  of  non¬ 
linear  svstems.  Generally  the  smaller  the  increment 
selected  for  h.  the  greater  the  accuracy  obtained 
in  the  res|»onsc  solution.  Kx|>eriencc  has  shown 
that  angular  inc  rements  0  =  360  cuA/2tt  between 
15  and  30  degrees  lor  si ngle-degree-of -freedom 
svstems  are  generally  acceptable  while  for  two- 
degree -of -freedom  systems  0's  should  range  be¬ 
tween  10  and  20  degrees.  The  two  examples  in 
this  re|M>rt  used  these  ranges  ol  values  for  0 
and  the  reader  can  see  the  difference  In- tween  the 
analog  computer  res|Mtnse  and  the  numerical 
tesults. 

While  the  examples  used  herein  are  mec  hanical 
models  of  structures  subjected  to  foundation 
motion,  the  numerical  method  is  applicable  for 
Imding  the  transient  res|M>nse  of  a  set  of  non- 
lineai  nonautonomous  differential  equations  for 
olhei  t\  |H*s  of  physical  systems. 
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APPENDIX  A 

MEAN  VALUE  DETERMINATION  OF  A  FUNCTION 


There  are  two  common  graphical  methods  for 
finding  the  mean  value  of  the  function  shown 
in  Fig.  A I  for  the  increment  At  =  A.  The  first 
method  is  to  draw  a  line  in  an  arbitrary  direc¬ 
tion  by  eye  such  that  the  sum  of  the  enclosed  areas 
between  the  curve  and  the  straight  line  is  zero. 
The  mean  value  of  the  function  is  the  distance 
from  the  abscissa  axis  to  the  straight  line  at  the 
middle  of  the  interval.  For  example,  line  cd  in 
Fig.  A 1  is  such  a  line  and  the  point  /  on  the  line 
determines  the  mean  value  F. 

The  second  method  consists  of  drawing  the  line 
ab  which  connects  the  end  points  of  the  curve 
segment.  At  the  middle  of  the  segment  measure 
the  distance  between  the  curve  and  the  line,  that 
is,  eg.  Measure  off  2/5  of  this  distance  above  (or 
below)  point  g  to  establish  the  mean  value  at  point 
/.  This  second  method  is  accurate  if  the  curve  is 
any  part  of  a  quadratic  parabola. 


Kig.AI  -  Determination  of  the  mean  value  of  a  function  Fit) 
lor  the  im  rrmrnt  A 


APPENDIX  B 

NUMERICAL  INTEGRATION  EQUATIONS 
FOR  LINEARLY  DAMPED  SYSTEMS 


Numerical  integration  equations  are  presented 
for  the  differential  equation 


The  applied  force  is  represented  by  the  para¬ 
bolic  method,  so  that 


y  +  2  atuy  +  a #*y  =  — 


(Bl) 


F  =  F.  +  Sn 


where  — »  <  a  <  *,  to  >  0.  The  region  of  a  is 
broken  into  seven  cases  as  billows: 


For  foundation  motion  Eq.  (Bl)  becomes 

x  +  2awx  +  to**  =  —  z .  (B2) 


(iase  I: 

a  =  0 

Case  II: 

0  <  a  <  1 

Case-  III: 

a  =  1 

Case  IV. 

a  >  1 

Case-  V: 

—  1  <  a  <  0 

Case  VI: 

a  =  —  1 

Case  VII: 

a  <  —  1 . 

By  means  of  similarity  between  Eqs.  (Bl)  and 
(B2),  it  is  only  necessary  to  change  the  follow¬ 
ing  in  the  numerical  integration  equations  for 
applied  forces  to  obtain  those  for  foundation 
motion.  let 


y»  -  xm,  i>„  =  Um,  Fm  =  — m ie'n  , 

A »  _  _  z w  S»  _  _  Sj  S| _ i  _  _  S|-i 

k  w*'  k  «**  k  w* 

If  the  input  is  given  by  the  foundation  velocity, 
which  is  represented  by  the  |>arabolic  method. 
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the  derivative  of  i  gives  the  foundation  accelera¬ 
tion  as 

••  S.  SI  ,  (2 1  l\ 

h  2  U  h) 

where 

S«  in  1 1  in 
Sn  l  =  iNt|  —  2  in  “H  2n  t- 


K(|tiations  (Bit  and  (B2)  are  solved  lor  each  of 
the  seven  cases. 


CASE  I:  a  =  0 

Input  -  Applied  Force 
or  Foundation  Acceleration 


yn.  I  =  y»  COS  d  +  C,  sill  ®  H - ^  (I  —  COS  d) 


CASE  II:  0  <  a  <  I  (Let  r  =  VI 

Input  -  Applied  Force  or 
Foundation  Acceleration 

yn.  i  =  y,  r  °®  |cos  rd  +  ~  sin  rflj 


+  p,  f  "® 


sin  rd 


+  I  —  **~"®  |(os  rd  +  *  sin  rdj  j 


+  4s  1  —  <  1  —  e  “®  cos  rd) 

K  U 


(  1  —  2al)  t~ 


^sitw^l 

rd  J 


si  , 

4a 

r 2(  1  -  4a* ) 

.2«1 

2* 

0 

l  e* 

e  J 

( 1  -  e'°*  cos  rd) 


1  -  2a»  2a (3  - 
d  d * 


~  4a1 ) 1 

•*  J 


(\  sin 

,  M  ■ 

[sin  d 

2(1—  cos  d 

r  ».  sin  r®] 

l1  d  ) 

[  » 

d* 

'  r  1 

=  — y,  sin  d  +  t>»  cos  0 

.  A  n  •  „  ,  Sn  / 1  ~  COS  d\ 

+  T  *,n  •  +  T  l  0 — ) 

S|  i  /I  +  cos  (J  2  sin  d\ 

2k  v  »  tj«  r 

Input  -  Foundation  Velocity 

. . „  .  ..  ■  „  s,  n  -  cos  d\ 


*n  COS  d  + 


„  5,  / 1  —  cos  0\ 

Un  Sltl  d - - - - 

CU  \  d  f 


.  sin  rd  J  „  a  .  \ 

•’»  •  i  =  — V  n  e  °® — - - He,  e  “®f  cos  rd-— sin  rd  I 

+  r  a.  *'» 

I  r 

S,  / 1  ~  f  °*  cos  rd  _  ar  a*  sin  rd\ 
k  \  d  rd  ) 

+s:k{!-(T+')rr7"’") 

f 2(1—  2a*)  al  r  "®  sin  rdl 

l  •*  r  J’ 

Input  •  Foundation  Velocity 


Si  ,  / 

1  +  c  os  d 

sin  d\ 

=  “•(cos  r0  +  —  sin  rtf) 

2d 

d1  / 

\  r  / 

S,  /l  —  e  “*  cos  rd  a  r*®  sin  rd 


S,  sin  d 


^_sin_rd\ 
rd  / 


n  .  „  on  .sill  i 

Un  *  I  —  —  *«  Sin  d  +  Un  COS  d - — 

cud 

_  Sn  I  /  I  ~  COS  d  _  Sill  d\ 
~  t o  \  d*  2d  / 


SI  ,  [1  /I  ,  2a\/l  -r  “•  cos  rd\ 


.  fa  (I  —  2 a* )  1  e  ***  sin  rd 
12  d  I  rd 
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...  sin  r9  ,  „  a  . 

U„.l  —  Xu  €  of - - - h  <>  “'IcOS  r“~~  s,n  ) 

_  sin  rfl  Sj- 1  T  1  -  e~at  cos  rff 

<o  e  rO  u  [  9 * 

_/ a  ]\  f  °*  sin  rO  1 

\e  2}  re  y 


CASE  III:  a  =  1 

Input  -  Applied  Force  or 
Foundation  Acceleration 

y..i  =  y«(l  +  6)  e~*  +  v»9  e * 

+  j1  tl  -  (1  +  *)  e*] 


CASE  IV:  a  >  I  (Let  q  =  \V  -  1) 

Input  -  Applied  Force  or 
Foundation  Acceleration 

y n ,  i  =  y„  e  aB  ^cosh  qO  +  ^  sinh  qd'j 

.  sinh  qS 
+  v„  e  a* - *- 

q 

+  -~j7  j^l  —  t  “®  ^cosh  qO  +  —  sinh  gfljj 

+  T  [  1  ~  T  ( 1  "  e  “>sh  q6) 


r  2  i 

^Sl  , 

(  4a 

[2(1—  4a*)  2al 

[d  +e  #)  (1  -e-*)j 

+  ~2k~ 

l  o 

L  o*  « J 

k..i  =  -y»Bt~*  +  *'■(!-  9)e  *  +  Be  * 

Input  *  Foundation  Velocity 

a,.,  =x„(l  +  »)<•*+ u.flr* 

(HH 

9r  *  +  u„<l  -  9)' -•  r  * 

CO 


( 1  —  e  a*  cosh  q9) 

.  fl  -2o*  ,  2«(3-4«*)l 
+  l  9  +  9*  J 

Slllh  ?g| 


«’».  I  =  -y»  r 


sinh  qO 


+  v,  r'°*  ^cosh  q9  —  ^  sinh  q9 ^ 


+  Eze  *•  ?inh..g* 


+ 


U.  ,  I  =  Xu 


5*  /l  —  «*  “*  cosh  td  _  ae~a*  sinh^O\ 

*  \  »  qO  ) 


_  [2(1-20*)  «1 

L  #*  9\ 

f  -  sm.h  ?.gj 
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Input  •  Foundation  Velocity 

r„,i  =  x„  f  "*  ^rosli  qtt  I  ^  sinh 
sinh 


+  u„  e 


S„  ^  1  —  e  “*  <  osh  qtt  _  ae  °*  sinh  q6 


'  sinh  q0\ 
tt  qO  ) 


S|_,  |  I 

h) 


U 

.2 


(1  ~  2«») 
tt 

sinh  qtt] 

qO  I 


CASE  VI  (continued) 

+f  [> +«•+!('-*')] 


Vn  ♦  i  = 


-y,0e»  +  v,(l  +  9)e* 

♦*-  +  W-) 


U  n  *  l  Xu  t3 


sinh  qtt 


■+  «»  **  "*  ^tosli  qtt  —  —  sinh 


<u  qtt 


-(H) 


Midi  gtfj 

<?»  J 


CASE  V:  -I  <  a  <  0 


Input  -  Foundation  Velocity 

x,,i  =  *»(  1  -  0)e*  +  Untte* 


SI  , 

1  -  r  "*  tosh  qtt 

SI  .  [1 

('  2\( 

l-e*\ 

0) 

tt1 

at  (.« 

\2  9/[ 

9  ) 

-(H)"] 


S. 


Use  lilt*  equations  lor  (;ise  II  but  replace  a  by  u„,  t  =  —xm9r'  +  u„(  1  +  0)e*  —  —  e# 

-a. 


CASE  VI:  a  =  -I 

Input  •  Applied  Force  or 
Foundation  Acceleration 

y. .  i  =  y»  ( I  -  9)r*  +  v*9r' 

+  Y  [I  -  (I  -#)<■"] 


CASE  VII:  a  <  -1 

Use  I  he  equations  lor  C.ase  IV  but  replace  a  by 

—a. 


